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Abstract 

During a core-collapse supernova, absorption of v e emitted from the proto-neutron star by pro- 
tons in the hydrogen envelope produces neutrons and positrons. Neutron capture on protons and 
positron annihilation then produce 7 rays of 2.22 and 0.511 MeV, respectively. We calculate the 
fluxes of these 7 rays expected from a supernova with an 11 Mq progenitor. The flux from neutron 
capture on protons exponentially decays on a timescale of 564 s, which is determined by neutron 
decay and capture on protons and 3 He nuclei. The peak flux is 2.38 x 10 -7 cm -2 s _1 for a super- 
nova at a distance of 1 kpc. In contrast, the 7-ray flux from positron annihilation follows the time 
evolution of the v e luminosity and lasts for ~ 10 s. The peak flux in this case is 6.8 x 10~ 5 cm -2 s _1 
for a supernova at a distance of 1 kpc. Detection of the above 7-ray fluxes is beyond the capability 
of current instruments, and perhaps even those planned for the near future. However, if such fluxes 
can be detected, they not only constitute a new kind of signals that occur during the gap of several 
hours between the neutrino signals and the optical display of a supernova, but may also provide a 
useful probe of the conditions in the surface layers of the supernova progenitor. 
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I. INTRODUCTION 



On exhaustion of nuclear fuels, the core of a massive star (> 8 M® with M® being the 
mass of the sun) undergoes gravitational collapse, thereby initiating the supernova process. 
Two classes of signals are expected from such an event: the neutrinos emitted by the proto- 
neutron star formed from the collapsed core and the photons radiated as the supernova shock 
emerges from the stellar surface. The neutrino signals start immediately after the shock is 
launched and last for ~ 10 s. However, it takes several hours for the shock to emerge, and 
consequently, the associated photon radiation, in particular the optical display, is delayed 
from the neutrino burst by this shock propagation time. In this paper we consider a third 
class of signals that occur before the shock emergence. These signals are 7 rays induced by 
neutrino reactions in the stellar envelope. 

We assume that the star undergoing core collapse still has its hydrogen envelope. As the 
neutrinos from the proto-neutron star stream through this envelope, the reaction 

v e + p — > n + e + (1) 

produces a neutron and a positron. Subsequently, the capture of the neutron through 

n + p-^ 2 H + 7 (2) 

produces a 7 ray of 2.22 MeV, and the annihilation of the positron can produce two 7 
rays of 0.511 MeV each. The above mechanism of 7-ray emission from supernovae has 



been considered earlier by Refs. 



2|. However, these studies only estimated the expected 
7-ray fluxes without giving an analysis of all the physical processes involved in the 7-ray 
production. For example, the thermalization of neutrons and positrons in the stellar envelope 
was not discussed, and neither was the detailed time structure of the 7-ray emission. We note 
that 7 rays from neutron capture on protons and positron annihilation were also discussed 
in the context of solar flares (see e.g., Ref. [sj]) and interstellar medium (see e.g., Ref. Q]), 
and there were extensive studies of the physical processes involved in the 7-ray emission 
from positron annihilation (see e.g., Ref. [5J). 

We here present detailed analyses of the major physical processes that lead to 7-ray 
production following the reaction v e + p — > n + e + in the hydrogen envelope of a massive 
star. In particular, we show that the 7-ray emission due to neutron capture on protons 
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lasts for ~ 10 3 s while that due to positron annihilation follows the time evolution of the 
v e luminosity and lasts for ~ 10 s. For concreteness, we adopt a specific model of neutrino 
emission and a specific stellar model for the conditions in the hydrogen envelope. Our 
analyses can be easily generalized to other neutrino emission and stellar models. 

We assume that the gravitational binding energy Eb of the final neutron star is emitted 
equally in u e , u e , v^, z/ M , v T , and v T and that the v e luminosity L Pe (t) decays exponentially 
on a timescale r. Thus, 

M()3(«pH/r). (3) 

We take Eb = 3 x 10 53 erg and r = 3 s. The normalized v e energy spectrum is taken to be 

1 E- 

UEpJ = ^F 2 (^ e )ex P [(^ e /T p J-^ e ] + l' (4) 

where T De =3.76 MeV, r] Pe = 3, and 

Fn{v) = / 7 r—rdx. (5) 

J exp(x - rj) + 1 

The corresponding average v e energy is (E 9e ) = Tp e F 3 (r] Pe ) / F 2 (?7p e ) = 15 MeV. The cross 
section for the reaction v e + p — > n + e + (see Appendix |A]) averaged over the above spectrum 
is ( a _ ep ) = 1.87 x 10" 41 cm 2 . 

For the conditions of the hydrogen envelope, we use the model of an 11 M star in 
Ref. [(J. The region of interest for 7-ray emission is limited by the interaction of 7 rays 
with matter. For 7 rays of ~ 1 MeV produced in the hydrogen envelope, the dominant 
interaction is Compton scattering on electrons. The relevant cross sections (see Appendix |A|) 
are cr 7 ( np)e = 1.38 x 10~ 25 cm 2 and a 7 ( e ±) e = 2.87xl0~ 25 cm 2 for 7 rays of 2.22 and 0.511 MeV, 
respectively. The corresponding mean free path is 

1 -1 aa , n 9 ( 10_8 S cm ~^ ( iO -25 cm 2 



1^ = = 1-66 x 10 9 % cm, (6) 

n e a-y e V pY e J \ a 7e J 

where n e = pY e NA is the electron number density, p is the matter density, Y e is the number 
of electrons per nucleon, and is Avogadro's number. The surface zone of our adopted 
stellar model has p = 1.59 x 10 -8 g cm -3 and Y e = 0.85, for which Z 7 ( np ) e = 8.9 x 10 8 cm and 
^7(e±)e — 4.28 x 10 8 cm for 7 rays of 2.22 and 0.511 MeV, respectively. The radius of this 
zone is R = 2.36 x 10 13 cm. For considering the emission of 2.22 and 0.511 MeV 7 rays (we 
do not treat the emission of scattered 7 rays at other energies), it is sufficient to focus on the 
outermost region with a radial thickness d satisfying l ie <^ d R. The stellar conditions 



stay constant in this region. In addition, the rate for production of neutrons and positrons 
per nucleon by the reaction v e + p — > n + e + is the same throughout this region and is 



where Y p = 0.7 is the number of protons per nucleon in the region, and Ap eP (0) = 1.30 x 
10~ 12 s _1 for the adopted parameters. 

We study the emission of 2.22 MeV 7 rays from neutron capture on protons in Sec. [Ill 
and that of 0.511 MeV 7 rays from positron annihilation in Sec. IIII1 We discuss our results 
and give conclusions in Sec. [TV] 

II. 7- RAY EMISSION FROM NEUTRON CAPTURE ON PROTONS 

Before discussing the physical processes leading to 7-ray emission from neutron capture 
on protons, we give a simple estimate of the time evolution of the corresponding flux. This 
evolution is closely related to that of the neutron number per nucleon Y n (t) in the stellar 
surface region from which 7 rays can escape efficiently The increase in Y n (t) is due to v e 
absorption by protons with the rate Xp eP {t) given in Eq. (JTj). The decrease in Y n (t) is caused 
by neutron decay and capture onto protons and 3 He. The neutron lifetime is 77 = 887 s. For 
capture of low-energy neutrons, the cross section is inversely proportional to the neutron 
velocity v n . Consequently, the product of the cross section and v n is independent of the 
neutron velocity distribution. We use (v n a np ) = 7.32 x 10~ 20 cm 3 s _1 for capture onto 
protons and (v n a n s) = 1.17 x 10~ 15 cm 3 s _1 for capture onto 3 He. The corresponding 
capture timescales are r np = {pY p N A^nanp))^ 1 = 2.04 x 10 3 s and r n3 = (pY 3 N A^n^nz)) 1 = 
6.38 x 10 3 s, where we have taken the number of 3 He per nucleon to be = 1.4 x 10~ 5 , the 

n 

same as estimated for the solar photosphere |3J . Note that the dominant channel of neutron 
capture onto 3 He is n + 3 He — >• p + 3 H, which does not produce any 7 ray. However, the 
enormous cross section of this channel enables it to play a significant role in determining 
the evolution of Y n (t) in spite of the small abundance of 3 He. 
A simple estimate of Y n (t) can be obtained from 





dY. 



n 



Y n (t) 



(8) 



dt 



4 



where r e g = r n 1 + r np l + r n ^ and r c fj = 564 s. The solution to the above equation is 

Y n (t) = A PeP (0)r (^—) [exp(-t/r eff ) - exp(-*/r)]. (9) 

As r cS ^> r, Y n (t) rises to its peak value 

Yn k ~ W0)r = 3.90 x 10~ 12 (10) 

on a timescale of ~ 10 s and then exponentially decays on the timescale r e ff. 

Due to Compton scattering, only those 7-rays of 2.22 MeV produced in the outermost 
stellar layer with a thickness of ~ / 7 ( np ) e — 8-9 x 10 s cm will escape efficiently The corre- 
sponding flux at a radius r > R can be estimated as 

-r , . s pY n (t)N A R 2 l^ np)e 7 /lkpc\ 2 _ 2 _ x 

$ 7 („ p )(r,* r ) — 2 ~ 9.53 x 10 exp(-t/r cff ) cm s \ (11) 

Tnp r \ f J 

where t r is the time at which the 7 rays emitted at time t arrive at radius r and we have 
used Y n (t) ~ y^ k exp(— t/r cS ) in the second approximation. To show the dependence on the 
model of neutrino emission and stellar conditions, we rewrite the above equation as 

•^^^-ssfefS)^-)- (12) 

It can be seen that the exact form of L 9e {t) is unimportant so long as the timescale of 
neutrino emission is r c g. In addition, the density of the stellar surface region controls the 
peak magnitude and the decay timescale of the flux via r np and T e g , respectively. Integrating 
the flux over time, we estimate the total fruence of 2.22 MeV 7 rays at radius r as 

~ ^iE~) [y^:) {-J ~ M? X 10 (—J Cm ' (13) 
Note that Fyfap) i n general still depends on the density of the stellar surface region due to 
the competition between neutron decay and capture. This dependence ceases only when 
neutron decay can be ignored (i.e., r n ^> r np ). 

A. Thermalization and Diffusion of Neutrons 

As discussed above, the timescales for neutron capture onto protons and 3 He are r np = 
2.04 x 10 3 s and r n3 = 6.38 x 10 3 s for the adopted density and composition of the stellar 
surface region. We now show that these timescales and the neutron lifetime are so long 
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that neutrons are thermalized due to scattering by protons before being captured or decay. 
Following the absorption of v e by protons, positrons are emitted approximately isotropically 
with energies of ~ E De . Using the v e energy spectrum in Eq. (J4J) and the cross section 
in Eq. flAip . we obtain the average recoil energy of the neutrons produced along with the 
positrons as (E™ c ) ~ [F 6 (77^ J /F 4 (?7y J ]T| e /M n ~ 543 keV, where M n is the neutron rest 
mass. The region of interest has a temperature T = 1.53 x 10 4 K corresponding to a 
thermal energy E t h = (3/2)fcT = 1.98 eV <C (E^ c ) with k being Boltzmann's constant. 
Consequently, neutrons lose energy through scattering by thermal protons until neutrons 
are thermalized (scattering by other particles can be ignored). The average decrease in the 

n 

natural logarithm of neutron energy per scattering is unity [7|: (\n(E n j + i/ E n j)) = — 1, 
where E n j is the neutron energy after j scatterings. So we approximately have 

E nd ~(E™)eM-j)- (14) 

The mean free path between scattering is 

1 7 /1(T 8 g cm~ 3 \ AO b\ . . 

he = = 1-66 x 10 7 %. cm, 15 

n p cr sc \ pY p J \ a sc J 

where n p = pYpN^ is the proton number density and <r sc is the scattering cross section. For 

the relevant neutron energies, a sc ~ 10 b [8| corresponds to / sc ~ 1.49 x 10 7 cm for the 

adopted stellar conditions. The timescale for thermalization can then be estimated as 

Jmax 1 ^ 

^therm ^ > ; ~ 15 S < (16) 

j=0 \/2>E ntj IM n 

where j max ~ H(E r n cc ) / E th ) ~ 13. 

Once thermalized, neutrons diffuse until they decay or are captured or escape from the 
stellar surface. The mean speed of thermal neutrons is v n = ^/8kT/ {jiM n ) = 1.80 x 
10 6 cm s^ 1 . If they escape from the star under consideration (11 M Q with a radius R = 
2.36 x 10 13 cm), they will not fall back onto the star as they will decay while they are still 
moving away from the star. Therefore, those neutrons that can escape will not contribute 
to the production of 7 rays. In order to escape, neutrons must diffuse to the stellar surface 
on timescales shorter than r e g . The thickness 5 of the layer from which neutrons can diffuse 
to escape can be estimated from 5 2 ~ (w n r e fr// SCit h)/s C t h> which gives 5 ~ ^/v n T cS l sc th = 
8.6 x 10 7 cm with i S c,th = 7.28 x 10 6 cm corresponding to cr SCit h = 20.5 b for thermal 
neutrons. As 5 ~ 0.1l 7 ( np ) e , the reduction of potential 7-ray production due to the escape 
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TABLE I: Rate, mean free paths, and timescales for the important processes involved in the 7-ray 
emission from neutron capture on protons. The adopted conditions in the stellar surface region 



are characterized by R = 2.36 x 10 13 cm, p = 1.59 x 10 8 g cm 3 , Y e = 0.85, Y p = 0.7, and 
y 3 = 1.4 x 10~ 5 . 



Ap eP (0) 


^y(np)e 


^sc,th 


Trip 


T n 3 


Ted 


1.30 x 10~ 12 s" 1 


8.9 x 10 8 cm 


7.28 x 10 6 cm 


2.04 x 10 3 s 


6.38 x 10 3 s 


564 s 



of neutrons is insignificant. For the outmost layer of thickness Z 7 ( np ) e from which 7 rays can 
escape efficiently, essentially all the neutrons in this layer decay or are captured by protons 
or 3 He during diffusion. Indeed, the timescale for neutrons to diffuse out of this layer is 

r diff = = 6.04 x 10 4 s » r cS . (17) 

For the adopted stellar conditions, r n = r np /2.3 and r n3 = 3.1r np . So a fraction r e fj/r np = 
[1 + (r„p/r n ) + {jnp/Tn^Y 1 = 27.6% of the neutrons are captured by protons to produce 
2.22 MeV 7 rays. The reduction of potential 7-ray production is mainly due to neutron 
decay and to a smaller extent due to neutron capture by 3 He. The rates and mean free 
paths for the important processes involved in the 7-ray emission from neutron capture on 
protons in the star under consideration are summarized in Table HI 



B. Evolution of Neutron Density and Emergent 7-Ray Flux 

Based on the preceding discussion, we consider the following quantitative model for cal- 
culating the flux of 2.22 MeV 7 rays from the star under consideration. We focus on the 
outermost layer of thickness d = 10 10 cm ^> l 7 r np ) e to find the emergent flux <J> 7 ( np ) (i?, t) at 
the stellar surface. As d R, we can treat the layer of interest as a slab perpendicular 
to the x-axis with the stellar surface at x = d (see Fig. [1]). Knowing the neutron number 
per nucleon Y n (x, t) in this layer, we can calculate the local flux of 2.22 MeV 7 rays in the 
positive x-direction as 

47rr np J J Q Ar np 
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FIG. 1: Sketch of the stellar surface region of interest to 7-ray emission. As L < d < R, this 
region can be treated as a slab of thickness d perpendicular to the x-axis with the surface at x = d. 
The calculation of the local 7-ray flux in the positive x-direction involves integration over the 
forward solid angle defined byO<0<7r/2. 

where 9 and <fi are respectively, the polar and azimuthal angles with respect to the x-axis. 
Consequently, the emergent flux at the stellar surface is 

<S>~ ( (n P )(R,t) = j^- [ F n (x,t)exp (-j — -\ dx, (19) 

which is related to the flux at r > R as $ 7 ( np )(r, t T ) = (i?/r) 2 $ 7 ( np )(_R, t). 
The evolution of Y n (x,t) is governed by diffusion: 

dY n d 2 Y n YJx,t) , . 

at ox 1 r cff 

where D = /scth^n/3 is the diffusion coefficient of thermal neutrons. As T e g is much 
longer than the timescales for neutrino emission and thermalization of neutrons, we as- 
sume Y n (x, 0) = Y£ k [see Eq. ([TO]) ]. Further, on the timescale relevant for 7-ray production, 
diffusion has little effect on the spatial distribution of neutrons at x = [see Eq. (TTTT) ] . So we 
take (dY n /dx) x= Q = as the inner boundary condition. For the outer boundary at x — d, we 
assume that the drift neutron flux given by diffusion is the same as the escaping neutron flux: 
-D(dY n /dx) x=d = v n Y n (d,t)/4, which is equivalent to l sc ,th{dY n /dx) x=d = -(3/4)Y n (d,t). 
With the initial and boundary conditions discussed above, the solution to Eq. (1201) is 

Y n (x, t) = 2Y? £ , Si *M) C ° s( "f exp[-(«jD + r^t], (21) 
Kjd + sm[Kjd) cos(Kjd) J 

where Kj satisfies Kjl sc ,th tan(Kjd) = 3/4 as required by the outer boundary condition. The 
evolution of Y n (x, t) given by Eq. (I2T1) is shown in Fig. Ett. 

Substituting the above expression of Y n (x,t) in Eq. ( |T9l) . we obtain 



*tW(R, t) = ^^nrte A . exp [-(^ + T £)t] 

& 1 rt.n 



(22) 



n P j= o 



where 

_ sin(/t 3 -rf)[cos(Kjd) + (ftj/ 7 (n P )e) sin(Kjri) - exp(-d/f 7 ( np)e )] 
J + sin(« J -d) cos(/Cjd)][l + (Kj/ 7 ( np ) e ) 2 ] 

As i S c,tii Z 7 (np)e *C it can be shown that 

A i ~ n , iA.,1 /o^in , z NaT- </ 4 J 



:,thJ 



The outer boundary condition gives jir/d < Kj < [(j + (l/2)]ir/d. For j <C 3d/(47rZ SC) t 
328, Kj ~ [j + (l/2)]7r/d. This approximation can be used for all j as there is little difference 
between j and j + (1/2) for sufficiently large j. Noting that 4«jZ S c,th/3 ~ [j + (l/2)]/328, 
Ki^rMe « b" + (l/2)]/3.58, and K |Dr efI « {[j + (l/2)]/64.1} 2 , we have 

*7(np)(n*r) = ( — j $ 7 (np)(-R,t) 

~ Y n k P N AR%{np)e , ,/ ^(np)e/d 

2r np r 2 6XP( VTeffj ^l + {b- + (l/2)M 7(np)e /dP (25) 

« 2.38 x 1CT 7 (^^— ^) exp(-t/r cfT ) cm" 2 s -1 . (26) 

The result in Eq. (1261) is obtained by approximating the sum in Eq. (1251) by the integral 
J* °°(l + y 2 )~ 1 dy/n = 1/2. This result is smaller than the estimate in Eq. (fTTj) by a factor of 
4, which comes from the integration over the solid angle to obtain the local flux at a point in 
the region of 7-ray production [see Eq. ( fl8l) ]. We numerically evaluate the flux of 2.22 MeV 
7 rays at r = 1 kpc as a function of time t from Eqs. ([221 and ([23]) and show the result in 
Fig. [2b along with the approximation in Eq. (T261) . The numerical result is essentially the 
same as the approximation. 



III. 7-RAY EMISSION FROM POSITRON ANNIHILATION 

The processes leading to the emission of 0.511 MeV 7 rays from annihilation of positrons 
subsequent to their production by u e absorption on protons are much more complicated than 
those resulting in the emission of 2.22 MeV 7 rays from neutron capture on protons. This is 
because a positron can lose energy through many processes such as electronic excitations of 
ionized and atomic matter, bremsstrahlung, Compton scattering, and synchrotron radiation 
in the presence of a magnetic field. In addition, a positron can directly annihilate with or 
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FIG. 2: (a) Neutron number per nucleon Y n (x,t) as functions of x in the stellar surface region for 
t = 500, 1000, and 2000 s. Note that the decline of Y n (x,t) with x due to diffusion only occurs 
in a narrow region near the surface, (b) Time evolution for the expected flux of 2.22 MeV 7 rays 
from neutron capture on protons in a supernova at a distance of 1 kpc as calculated by solving the 
diffusion equation (solid curve) and estimated by neglecting diffusion (dashed curve). 



form a positronium (Ps) with a free or bound electron: 



e + + e" 



e + + H 
e + + He 
e + + e~ 

e + + H 
e+ + He 



7 + 7, 

H+ + 7 + 7, 

He + + 7 + 7, 

Ps + 7, 
Ps + H+, 
Ps + He+. 



(27) 
(28) 
(29) 
(30) 
(31) 
(32) 



In the case of Ps formation (with free and bound electrons through radiative combination and 
charge exchange, respectively), the singlet Ps ( x Ps with total spin 0) decays by emitting two 
photons of 0.511 MeV each while the triplet Ps ( 3 Ps with total spin 1) decays by emitting 
three photons with a continuous energy spectrum. Here we focus on the production of 
0.511 MeV 7 rays only. We will loosely refer to both direct annihilation of positrons and Ps 
decay as "annihilation" of positrons. 

The energy loss and the annihilation processes of the positrons depend on the ionization 
state of the medium. We consider that the stellar surface region is in local thermodynamic 
equilibrium at T = 1.53 x 10 4 K and p = 1.59 x 10~ 8 g cm" 3 . We assume that the forward 
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and reverse processes of the following chemical reactions are in equilibrium: 

e- + H+ ^ H + 7, (33) 
e" + He+ ^ He + 7. (34) 

As we will see shortly, nearly all of the He atoms remain neutral and therefore, we can 
ignore the presence of He ++ . To avoid confusion, we denote the numbers of free e~ , H + (free 
protons), He + , and neutral H and He atoms per nucleon as Y e -, Y H +, Y ne +, Yr, and Y Rei 
respectively. These quantities satisfy 

F H + +Y K = Y p , (35) 
Y Re+ + F H c = Y a , (36) 
Y H+ + Y He+ = Y e -, (37) 

where Y p = 0.7 and Y a = 0.075 are the numbers of protons and 4 He nuclei per nucleon, 
respectively (Y p + 4Y a = 1 for consistency with a medium where the dominant nuclei are 
protons and a particles). 

Based on the chemical equilibrium for the reactions in Eqs. ( 1331) and (1341) . we obtain 



(38) 



Y He+ 1 f2g Re+ \ fm e kT\ 3/2 f _ 1.52 x 10~ 2 

Y Uc ~ P Y e -N A Uhc J V 2vr^ J 6XP \ kT J ~ Y e - ' [6J) 

where m e is the electron rest mass, gn ~ 4.8, g He + ~ 2, and gn e ~ 1 are the partition 
functions of the corresponding species, and /h = 13.6 eV and I Re = 24.6 eV are the first 
ionization potentials for the ground states of the H and He atoms, respectively. There is 
some subtlety in calculating the partition functions as discussed in Appendix El However, 
this does not affect the general results: nearly all of the H atoms are ionized and nearly all 
of the He atoms remain neutral. Solving Eqs. (!35|) - (l39|) . we obtain Y e - 0.67, Y a + rs 0.667, 
Y R w 0.033, Y Rc w 7.33 x 10" 2 , and Y He+ w 1.7 x lO" 3 . 

The positrons produced by v e absorption on protons have an average energy (E e +) ~ 
[^(r/pj / F i (rjp e )]T Pe ~ 21.2 MeV. As they pass through the essentially ionized plasma of 
the stellar surface region, they can lose energy through many processes, directly annihilate, 
and form Ps. The x Ps formed decays into two photons with an extremely short lifetime of 
t 27 = 1.25 x 10~ 10 s while the 3 Ps decays into three photons with a much longer lifetime of 
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T3 7 = 1.42 x 10 7 s. Consequently, once formed, 1 Ps immediately decays while 3 Ps can be 
broken up or converted into 1 Ps by the following reactions before it decays: 



e~ + 3 Ps -> e~ + e + + e", (40) 

H + 3 Ps -> H + e + + e", (41) 

H + + Ps -> H + e + , (42) 

He + + Ps -> He + e + , (43) 

e~ + 3 Ps -> e" + ^s, (44) 

H + 3 Ps -> H + ^s. (45) 



Note that the breakup reactions in Eqs. ( 1421) and ( 143|) are the reverse processes of the Ps 
formation reactions in Eqs. (j3Tj) and (1321 . Although all the processes in Eqs. fl27|) - fl32|) and 
(l40|) - (145l) are involved (to varying degrees) in the 7-ray emission from positron annihilation 
(see Fig. [3^), for the adopted stellar conditions, the net result is rather simple: (1) Before 
thermalization, w 12% of the positrons directly annihilate and « 4% form Ps. Based 
on the available spin states, 1/4 of the Ps formed are 1 Ps, which immediately decay into 
two 0.511 MeV 7 rays with line widths of ~ 6 keV, and 3/4 of the Ps formed are 3 Ps, 
which are immediately broken up (see Fig. [3b). The positrons released by breakup of 3 Ps 
are quickly thermalized. (2) So effectively ~ 87% of the initial positrons are thermalized. 
Essentially all of the thermal positrons form Ps. The 3 Ps formed are immediately broken 
up or converted into 1 Ps (see Fig. |3b). The positrons released by breakup of 3 Ps are again 
quickly thermalized. All the x Ps formed after thermalization immediately decay into two 
0.511 MeV 7 rays with line widths of ~ 2 keV. 



A. Thermalization of Positrons 



As mentioned above, positrons can lose energy through many processes. For the adopted 
stellar conditions, the dominant energy loss is due to excitation of the free electrons in the 
plasma. The thermalization of positrons spans the highly relativistic and nonrelativistic 
regimes. We use the general results on positron interaction with plasma electrons given in 
Ref. [9(. In the relativistic regime, the energy loss per unit length of propagation is 
'dE P 



dx 



ex, pi 



AnpY e -N A (-^] B rcl = 4.88 x 10~ 9 (-V y e -£ rel MeV cm" 1 , (46) 
\m e v z J \v/ 
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FIG. 3: (a) Summary of the processes involved in the 7-ray emission from annihilation of positrons 
with an initial energy of ~ 20 MeV. Direct annihilation with electrons before thermalization pro- 
duces 7 rays (indicated by 7*) with energies shifted far from 0.511 MeV. Positronium (Ps) formation 
with free and bound electrons is referred to as radiative combination and charge exchange, respec- 
tively, (b) Summary of the consequences of Ps formation that occurs after the positron kinetic 
energy drops below ~ 250 eV but before thermalization. (c) Summary of the consequences of Ps 
formation after thermalization. 



where 



B rci = In 



a/25(7 - 1) m e vc 



Tiw ri 



+ 6(7,*). 



(47) 



In the above equations, v is the positron velocity, c is the speed of light, 7=1/ a/1 — (v/c) 2 , 
5 is the maximum fraction of the positron energy lost in a single interaction and is taken to 
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be 1/2, oo p = a/ 4:7rpY e - Nac 2 / m e is the plasma frequency, and 



6(7,5) 



r 



i 



2 7 2 
7-1 



5 + 



l/ 7 -l 

7 



1 

2 V- 



1 



1 

+ 2 



7 



7 + 1 



7 + 2 

7 



1 1 3 

2 + 7 + V 



7 



7 Z 
d 2 



s2 + i 

1 



7 



7 



7 



7 

2 'a* 



4 
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The energy loss rate in the nonrelativistic regime is obtained by replacing B re \ in Eq. fT4T)|) 
with 



B m = In 



m P v 



e — 'Shpiiac/v) 



(49) 



where e ~ 0.577 is Euler's constant, a = e 2 /(hc) is the fine-structure constant, and ^iip{z) 
is the real part of the digamma function ip{z) = (dT/dz)/T(z) with T(z) being the Gamma 
function. 

The positrons produced by v e absorption on protons have an average energy (E e +) ~ 
21.2 MeV corresponding to (7) ~ 41.5. Once thermalized, they have an average kinetic 
energy E t ^ = (3/2)kT = 1.98 eV corresponding to 7th — 1 = 3.87x 10 -6 . Using the energy loss 
rate discussed above, we can estimate the total distance (Ax) th and time (At) th covered by a 
typical positron during thermalization. We proceed by noting that both — (dE e +/dx) eXiP \ and 
— (dE e +/dt) eXtP i = —v(dE e +/dx) e ^pi increase as the positron energy decreases. Consequently, 
the bulk of (Aa;) t h and (At) t h is covered before positrons become nonrelativistic, and we 
obtain 



(Ax] 



th 



dE P 



[E e+ ) 



3 x 10 8 cm, 



(50) 



{Ee+) (dE e +/dx) eXtP i 47rpY e -N A [e 4 /(m e c 2 )}B rel 

where we have taken advantage of the slowly-varying function B Te \ in evaluating the integral 
and used B Te \ ~ 20 for the adopted stellar conditions. The distance (Ax) t h estimated above 
corresponds to a thermalization timescale (A£) t h ~ (Ax) t h/c ~ 10~ 2 s. 

The energy loss rate discussed above is shown as a function of 7 — 1 in Fig. HI For 
comparison, the rates due to the other processes discussed in Appendix O are also shown. 
It can be seen that the energy loss due to excitation of the free electrons in the plasma 
dominates over the entire energy range between production and thermalization of positrons. 
We have numerically evaluated (Ax) t h and (At) th first using the actual energy loss rate due 
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FIG. 4: Rates of positron energy loss —dE e + /dx as functions of 7 — 1 for excitation of free electrons 
in the plasma (solid curve), excitation and ionization of bound electrons in atoms (dash-dot-dotted 
curve), bremsstrahlung (dashed curve), and Compton scattering on thermal photons (dashed- 
dotted curve). The rate for the first process is shown down to a positron kinetic energy of 2 eV 
(the thermal energy for the stellar conditions adopted here), that for the second process is shown 
down to a positron kinetic energy of 100 eV (the average excitation energy being 15 and 41.5 eV for 
H and He atoms, respectivley), and those for the last two processes are shown only for relativistic 
positrons with 7 > 3. Note that the actual rate multiplied by a factor of 10 3 is shown for Compton 
scattering. 

to excitation of the free electrons in the plasma only and then including the actual rates 
of the other processes discussed in Appendix The inclusion of the other processes has 
rather small effect and reduces (Ax)th from 2.98 x 10 8 to 2.49 x 10 8 cm and (A£) t h from 
10" 2 to 8.4 x irr 3 s. 



B. Direct Annihilation of Positrons before Thermalization 

In addition to losing energy, a positron can directly annihilate with a free or bound 
electron before being thermalized. As the thermalization distance (Ax)th is predominantly 
covered while the positron is still relativistic, direct annihilation before thermalization also 
predominantly occurs during the relativistic regime (see below). We ignore the 7 rays 
produced from such annihilation as the large positron velocity causes large shifts from the 
0.511 MeV line emitted in the center-of-mass frame for the annihilating positron and electron. 
The quantity of interest to us is the probability for direct annihilation of the positron before 
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thermalization. 

We first consider direct annihilation with free electrons. For a fast positron (v S> 2nac) 
annihilating with a free electron at rest, the cross section is 

The above result ignores the Coulomb interaction between the positron and the electron. 
This interaction becomes important at lower positron velocities, for which the cross section 

n 

[10J can be obtained by adding a "Coulomb focusing" factor to the low-velocity limit of the 
result in Eq. ( 15T1) : 

slow = ( e 2 V 2nMcM 2 ( , 2) 
da ' f \m e c 2 J 1 - exp(-2irac/vY 1 ' 

The probability for direct annihilation with free electrons before thermalization can be esti- 
mated as 

P<i*,t ~ Y e -pN A / — — / d-f ~ 0.1, (53) 

J {Ee+) {dE e +/dx) eXtP i 4:B rel J 1 7 

where we have considered only the relativistic regime and used the excitation of free elec- 
trons as the dominant energy loss mechanism with B ie \ ~ 20. Note that although for low 
positron velocities uf°J increases as ~ 1/f 2 , this is cancelled by the same increase of the 
corresponding —(dE e +/dx) eXiP \ [see Eq. (jlBI) ]. Thus, direct annihilation with free electrons 
before thermalization predominantly occurs when the predominant part of (Ax)th is covered, 
i.e., when the positron is still relativistic (see Sec. IIII Aj) . Note also that the probability of 
this occurrence depends logarithmically on the initial positron energy and the density of the 
medium (through B ie \). 

The cross section cr^fb f° r direct annihilation of a fast positron with a bound electron in 
an atom (or ion) is the same as <7^\ ■ In contrast, for v < ac, 

.2 \ 2 



a. 



slow _ _ rr I e \ C 
da,b — ^^eff I n , 



TiZ eS [ -, (54) 



where Z e g is a function of v and depends on the atom (see Appendix [D] for the case of the 
H atom) . The cross sections <jf°^i and o"daHe f° r direct annihilation with the electrons in 



the H [llj and He [12| atoms, respectively, are compared with erf™ in Fig. As <rf°J > 
"daH ' ff daHc an d K- ^ ^h, ^He, the probability for direct annihilation with bound electrons 
at low positron velocities is much smaller than the corresponding probability in the case of 
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FIG. 5: Cross sections for direct annihilation of slow positrons with free electrons (<7| .?Y , solid 



curve) atoms as functions of the positron velocity v in units of etc. 

free electrons, which is already very small. On the other hand, inclusion of the bound 
electrons increases the number of targets for direct annihilation in the relativistic regime. 
This gives P da « P dajf + P dajH + Pda,He ~ {Y e /Y e -)P da j, where P da is the total probability 
of direct annihilation before thermalization and Y e /Y e - = 1.27. Using the actual cross 
sections and rates of various energy loss processes over the entire thermalization regime, we 
have numerically calculated the probability of direct annihilation with free electrons before 
thermalization and found Pd a ,f = 9.52%. So the total probability including annihilation with 
bound electrons is P da ~ 12%. 

C. Ps Formation before Thermalization 

In addition to direct annihilation, positrons can form Ps with free and bound electrons 
before thermalization. As discussed above, direct annihilation predominantly occurs in 
the relativistic regime. As positrons become more and more nonrelativistic, Ps formation 
becomes more and more likely. The cross section for direct annihilation of slow positrons 
with free electrons multiplied by the number of free electrons per nucleon in the stellar 
surface region under consideration, Y e -af°J, is compared in Fig. [6] with the corresponding 
quantities F e -op Sj f, 1hO"p s ,h, and >HeCps,He for Ps formation with free electrons and the 
electrons in the H and He atoms, respectively. It can be seen that 1h<tp S! h and lH e cps,He are 
extremely large at 6.8 < E^f < 100 eV and 17.8 < E^f < 250 eV, respectively. As a result, 



curve) and with the electrons in the H (c d °H [Ul|, dash-dotted curve) and He (<r 
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da, He 



17 




1 1.5 2 

log(£j'7eV) 

FIG. 6: Cross section multiplied by the number of targets per nucleon for direction annihilation of 
slow positrons with free electrons (Y e -crf°J, s °hd curve) as a function of compared with the 
corresponding quantities for Ps formation with free electrons (Y e -a-p s j, dash-dottd curve) and with 
the electrons in the H (Ihcps.h, dash-dot-dotted curve) and He (I'HeO'Ps.He, dashed curve) atoms. 
The numbers of targets per nucleon used are Y e - = 0.67, Yh = 0.033, and Y^ e = 7.33 x 10~ 2 . The 
measured cross sections <7p Si H and <7p s hp for E^_ above the threshold of 6.8 and 17.8 eV but below 
100 and 250 eV are taken from Refs. |l3fl and respectively. The probabilities for Ps formation 
with the electrons in the H and He atoms are negligible for E + > 100 and 250 eV, respectively, 
before thermalization as <Jp s ,H and <Tp s ,He rapidly decrease for such E^. 

once positrons are slowed down to < 250 eV, they predominantly form Ps. Using the 
cross sections for Ps formation op S) f (see Appendix |D]), op Sj H |l3|, and op Si He U4 and the 
rates of various energy loss processes, we find that the probability for Ps formation with 
the electrons in the H and He atoms before thermalization is Pp Sj H ~ 1% and Pp s ,He ~ 3%, 
respectively, and that Ps formation with free electrons before thermalization is negligible 
(with a probability of -C 1%). So the total probability for Ps formation before thermalization 



is P Ps » P Ps ,h + Pp s . 



He 



4%. 



Based on the available spin states, x Ps and 3 Ps constitute 1/4 and 3/4 of the Ps formed, 
respectively. Thus, a fraction Pp s /4 pa 1% of the initial positrons form x Ps before thermal- 
ization and immediately decay into two 7 rays due to the extremely short lifetime of x Ps. 
The velocity distribution of the decaying 1 Ps gives rise to a spread in the energies of these 
7 rays, which are centered at 0.511 MeV. This spread can be characterized by a formal full 
width at half maximum (FWHM) of ~ 6 keV [5]. In contrast, the relatively long lifetime of 
3 Ps allows it to be broken up under the adopted stellar conditions. The dominant breakup 
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TABLE II: Probabilities and rates of various processes involved in the 7-ray emission from positron 
annihilation in the stellar surface region. The conditions in this region are characterized by p = 
1.59 x 1(T 8 g cm" 3 , T = 1.53 x 10 4 K, Y e - ps 0.67, Y R + ps 0.667, Y n ps 0.033, Y Rc ps 7.33 x 10~ 2 , 



and Y T 



He 4 



1.7 x 10 . Note that the rates for breakup and conversion of Ps depend on the 



formation mode of 3 Ps. The rates with the superscript "bt" correspond to the 3 Ps formed before 
thermalization and those with the superscript "Ps,H" correspond to the 3 Ps formed by the reaction 
e + + H — > Ps + H + after thermalization. 



Processes 



before thermalization 



after thermalization 



e + + e — ► 7 + 7 
e+ + H^H+ + 7 + 7 
e + + He -> He + + 7 + 7 
e + + e~ — > Ps + 7 
e+ + H -> Ps + H+ 
e+ + He -> Ps + He + 
e~ + 3 Ps — > e~ + e + + e" 
H + 3 Ps -> H + e+ + e~ 
H+ + Ps -> H + e+ 
He + + Ps -> He + e+ 
e" + 3 Ps -> e" + x Ps 
H + 3 Ps ->H + x Ps 



P da>f = 9.52% 
P da ,H « 0.47% 

P da)He « 2.08% 
Pp S ,f « 1% 
Pps.H ~ 1% 
Pp S ,He ~ 3% 

A^ 0)f » 4.5 x 10 8 s" 1 
A^ 0)H « 7.9 x 10 6 s" 1 

A 3^0.H+ ~ 7 ' 7 X 107 s_1 
^l ,He + « 9 - 8 >< 104 

A^ l f « 1.9 x 10 6 s" 1 

A 3^1.H ~ 3 - 9 X 106 s_1 



Ma,f 



7.7 x 10 2 s" 1 



Ada,H ~ 13 S 1 
Ada,He ~ 15 S _1 

A Ps f = 4.6 x 10 3 s" 1 
Ap s ,H ~ 9.9 x 10 4 s" 1 

Ap s ,He ~ 4.6 S" 1 



K% « 1.2 x 10? s- 1 
aEJh - 1.1 x 10 5 s" 1 

A 3^0,H+ ~ 3 - 2 X 108 S_1 



Af '" C H + ps 5.5 x 10 4 s- 1 
X P 3 tlf ~ 5.1 x 10 8 s- 1 
Ag^i H » 9.0 x 10 6 s" 1 



reactions are e~ + 3 Ps — > e~ + e + + e~ and H + + Ps — > H + e + with rates of ~ 4.5 x 10 8 
and 7.7 x 10 7 s _1 , respectivley [5||. The total rate for breakup of 3 Ps before thermalization 
is Agi^o « 5.3 x 10 s s _1 , to be compared with the rate of t^ 1 = 7.04 x 10 6 s _1 for 3 Ps 
decay. The rates of the above breakup reactions and many other reactions involved in the 
7-ray emission from positron annihilation are given in Table [IT] for the adopted stellar con- 
ditions. Thus, before thermalization a fraction (3/4)Pp s pa 3% of the initial positrons form 
3 Ps, which are immediately broken up. The positrons released by the breakup reactions are 
quickly thermalized. As a result, effectively a fraction 1 — P da ~ (-Pps/4) ~ 87% of the initial 
positrons are thermalized. 
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D. Fluxes of 0.511 MeV 7 Rays 



The thermal positrons again can directly annihilate with or form Ps with free or bound 
electrons. The rates for these processes are given in Table [IT] for the adopted stellar condi- 
tions. The dominant reaction is Ps formation with the electron in the H atom with a rate 
of Ap Si H = 9.9 x 10 4 s _1 [5]. The 3 Ps formed is quickly converted into 1 Ps or broken up. The 
dominant conversion reaction is e~ + 3 Ps — > e~ + 1 Ps with a rate of = 5.1 x 10 8 s _1 , 
and the dominant breakup reaction is H + + Ps — > H + e + with a rate of A 3 _^ — 3.2 x 10 8 s _1 
(5! (see also Table HT|) . Whether directly formed or produced by the conversion reaction, all 
ihe x Ps immediately decay into two 7 rays centered at 0.511 MeV with a FWHM of ~ 2 keV 
s|. The positrons released by the breakup reaction are again quickly thermalized and follow 
the same fate of thermal positrons as outlined above. 

Based on the discussion in the preceding subsections, a fraction Pda ~ 12% of the 
positrons produced by v e absorption on protons directly annihilate and a fraction Pp s « 4% 
form Ps before being thermalized, and the rest are thermalized on a timescale of ~ 10~ 2 s, 
which is much shorter than the timescale of r = 3 s governing the production of the initial 
positrons. The positrons released by the breakup of 3 Ps are also quickly thermalized. As the 
positrons are produced approximately isotropically, their spatial distribution can be taken 
as uniform before and after thermalization. So in the region of interest to 7-ray emission, 
the time evolution for the numbers of various species per nucleon is governed by 

^Vlp s bt L . ( , ^Ps.bt 

« -i 3 p s Ap eP (t) '—, (55) 

at 4 T2 7 

» 3 -P Ps X, eP (t) - A^ F3 Ps>bt , (56) 

dY e + ~ It td id \\ U\ , \bt \r , \Ps,H, 



dt 



,1 — -Pda — -Pp S )Ap e p(t) + A 3 ^ Y3 Ps!bt + A 3 ^ ^3Ps — Ap Sj HV e +, (57) 



dYi Ps 1 \Ps,h^ Yips , KQ s 

—7— ~ 7Ap B ,Hr e + + A 3 4 1 rsp B , (58) 

at 4 r 27 

~ ^A PSiH r e + - (A 3 ^ + A 3 ^ )Y3 Ps , (59) 

where Yip S) bt, isps,bt> ^Psj an d ^ 3 Ps correspond to the numbers per nucleon for the x Ps and 
3 Ps formed before and after thermalization, and Y e + is the number of thermal positrons per 
nucleon. As the timescale of r = 3 s governing Xp eP (t) is much longer than r 27 , (A 3 ^ ) _1 , 
(Aps,h)~\ and (A^ 1 , 1 + A 3 .!^ ) -1 , the above differential equations can be solved to good 
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approximation by setting all the time derivatives to zero. Thus, we have 

Yi PsM (t) « ^Pp s Ap eP (t)r 27 , (60) 
Yi Ps (t) ps M - P da - ip Ps ^) A PeP (t)r 27 . (61) 

These results are equivalent to what is stated at the beginning of Sec. IIIII and in Sec. IIII CI 
a fraction Pp s /4 ps 1% of the initial positrons form 1 Ps before thermalization and effectively 
a fraction 1 — P da — (Pp s /4) ps 87% of them form 1 Ps after thermalization. 

Using a similar prescription to that for estimating the flux due to neutron capture on 
protons, we find that at radius r > R the flux of 0.511 MeV 7 rays from decay of the 1 Ps 
formed before thermalization is 

y 1PsM (t) /pa 2 / d-aA, 

®7(e±),bt( r > *»•) ~ — : — / ex P -7 ax, (62) 

4r 27 V r / Jo V ^(e±)e/ 

3 J pN A l y{e ± )e , (63) 
exp(-t/r), (64) 



16 \r r 
P Ps E B ( Y p (a DeP ) 



3847rr 2 (P Pe )r \Y e a l{e ± )ei 
^ 7.8 x 10~ 7 0~^r\ ex P (-t/r) cm" 2 s" 1 , (65) 

where the factor 1/4 in Eq. ( 1621) comes from the integration of the local flux at a point in 
the stellar surface region over the forward solid angle. Note that although two 0.511 MeV 
7 rays are emitted in 1 Ps decay, only one contributes to the emergent flux (the other being 
emitted towards the stellar interior). Likewise, at radius r > R the flux of 0.511 MeV 7 rays 
from decay of the x Ps formed after thermalization is 

ps 6.8 x 10~ 5 ^i^^ exp(-t/r) cm" 2 s -1 . (67) 

The 7 rays associated with $ 7 ( e ±) i bt( r ; ^r) an d $ 7 ( e ±)(r, t r ) differ in that the former have a 
FWHM of ~ 6 keV while the latter have a FWHM of ~ 2 keV [5j. Note that $ 7(e ± ))bt (r, t r ) 
has only a weak dependence on the density of the stellar surface region through Pp s , which 
in turn depends logarithmically on the density through P nr associated with the energy loss 
due to excitation of free electrons [see similar density dependence for P da f through P rel as 
exhibited in Eq. (j5"3"]) ]. In contrast, apart from the weak density dependence of Pd a and Pp s , 
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$ T ( e ±)(r, t r ) is sensitive to the breakup of 3 Ps and the conversion of 3 Ps into 1 Ps, the rates 
of which are proportional to the density of the stellar surface region. For reasonable stellar 
conditions, $ 7 ( e ±),bt( r ? tr) is overwhelmed by $ 7 ( e ±\(r, t r ). 

IV. DISCUSSION AND CONCLUSIONS 

We have calculated the expected fluxes of 2.22 and 0.511 MeV 7 rays from neutron capture 
on protons and positron annihilation, respectively, following v e absorption on protons in the 
hydrogen envelope of an 11 M & star that undergoes core collapse to produce a supernova. 
The 7-ray flux from neutron capture on protons exponentially decays on a timescale of 
T c g = 564 s, which is determined by neutron decay and capture on protons and 3 He nuclei. 
The peak flux is 2.38 x 10~ 7 cm~ 2 s _1 for a supernova at a distance of 1 kpc. In contrast, 
the 7-ray flux from positron annihilation follows the time evolution of the v e luminosity. 
Although exponential decay on a timescale of r = 3 s is assumed here, the identical time 
evolution for the v e luminosity and the 7-ray flux from positron annihilation holds in general 
so long as the timescales for thermalization of positrons and Ps formation are much shorter 
than ~ 1 s. The peak flux in this case is 6.8 x 10~ 5 cm~ 2 s^ 1 for a supernova at a distance 
of 1 kpc. Detection of the 7-ray fluxes quoted above is beyond the capability of current 
instruments, and perhaps even those planned for the near future. For example, the proposed 
Advanced Compton Telescope jl5j has a spectral resolution of 0.2-1% over the energy range 
of 0.2-10 MeV and an angular resolution of ~ 1°, which are ideal for detecting the narrow 
7-ray lines discussed here. This instrument has a projected sensitivity of 5 x 10~ 7 cm -2 s _1 
for narrow lines but an exposure time of 10 6 s is needed. As the fluxes of 2.22 and 0.511 MeV 
7 rays discussed here only last for ~ 10 3 and 10 s, respectively, their detection requires much 
more sensitive instruments. 

If the 7-ray fluxes discussed here can be detected, they not only constitute a new kind 
of signals that occur during the gap of several hours between the neutrino signals and the 
optical display of a supernova, but may also provide a probe of the conditions in the surface 
layers of the supernova progenitor. For example, both the peak and the decay timescale 
r c fj of the 7-ray flux from neutron capture on protons depend on the density of the stellar 
surface region. A higher density decreases the timescales for neutron capture on protons 
( r np) and 3 He nuclei (173), which increases the peak flux and decreases r e ff. On the other 
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hand, a higher 3 He abundance decreases r n3 , and hence T e s, but does not affect the peak 
flux as neutron capture on 3 He nuclei consumes neutrons without producing any 7 ray. As 
another example, the rates for the breakup of 3 Ps and the conversion of 3 Ps into x Ps are 
proportional to the density of the stellar surface region. While the 7-ray flux from positron 
annihilation quoted above also applies approximately to higher densities than adopted here, 
for sufficiently smaller densities the 3 Ps formed would predominantly decay into three 7 rays 
with a continuous spectrum instead of being broken up or converted. This would decrease 
the flux of 0.511 MeV 7 rays from l Ps decay. In the limit where no 3 Ps are broken up 
or converted, this flux is reduced to 1/4 of the value estimated here. The above examples 
clearly illustrate that neutrino-induced 7-ray emission from the hydrogen envelope of a core- 
collapse supernova may serve as a useful probe of the conditions in the surface layers of the 
supernova progenitor. 
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APPENDIX A: CROSS SECTIONS FOR u e ABSORPTION ON PROTONS AND 
COMPTON SCATTERING 

In this appendix we set K — c — 1. The cross section for the reaction u e + p — > n + e + is 



Gp cos 2 6 C „ 2 



7T 



(f + 3g 2 )(l + 5 R )(E, e -Ay 1- 



2[p + 2(f + f 2 )g + 5g 2 } (E, 



9.56 x 10 



-11 / E i>e - A 

MeV 



, E 9 
1-7.2' 



M 



N 



cm 



M, 



N 



(Al) 



where Gf is the Fermi constant, 9c is the Cabbibo angle with cos^c = 0.9738, / = 1 and 
g = 1.27 are the vector and axial vector coupling constants, fi = 3.706 is the anomalous 
nucleon isovector magnetic moment, <5r ~ 0.024 is the inner radiative corrections, A = 
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M n — M p = 1.293 MeV is the difference between the neutron and proton masses M n and 
M p , and M N = (M n + M p )/2 = 938.9 MeV. 
The cross section for Compton scattering is 

2 \ 2 r 



7T 



e^y v 7; 2 e 7 2(1 + 2e^) 2 



(A2) 



where e is the magnitude of the electron charge, m e is the electron rest mass, and e 7 = i? 7 /m e 
is the photon energy in units of m e . 

APPENDIX B: PARTITION FUNCTIONS OF ATOMS 

The number of electron energy states in an isolated atom is infinite. This presents a 
problem in summing over these states to obtain the partition function of the atom as the 
sum formally diverges. However, the application of the partition function is sensible only 
when there are a large number of atoms. Consequently, the largest orbital radius of the 
electron in an atom is physically restricted to the interatomic distance. For our problem, 
the maximum radius r max can be estimated from 

A-jr 

T r 3 max p(Y p + Y a )N A =l. (Bl) 

Take the H atom as an example. The largest orbital radius of the electron is related to 
the maximum principal quantum number n max as r max ~ n max fi, 2 /(m e e 2 ). For our adopted 
stellar conditions, n max ~ 25. The partition function of the H atom is then 



g H = A^n 2 exp 

71=1 



kT \ n 2 



4.8, (B2) 



where the factor of 4 comes from the spin states for the proton and the electron, and the 
factor n 2 accounts for the degeneracy of the orbital states of the nth energy level. 

It can be seen that the partition function of the H atom under the stellar conditions of 
interest is dominated by the contribution from the ground state. This is also true for the 
He + ion and the He atom. As these two species are minor components in the stellar region 
of interest, we take their partition functions to be given approximately by the contributions 
from the corresponding ground states only, i.e., g He + ~ 2 and gn c 1. 
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APPENDIX C: RATES OF POSITRON ENERGY LOSS 



In addition to excitation of the free electrons in the plasma, positrons can also lose energy 
through excitation and ionization of the bound electrons in atoms and ions. Using the results 
from Ref . [l?| , we find that the energy loss rate per unit length of propagation for the latter 
process is 



dE e ^ 
dx 



AnpN A 



cx,at 



x 10" 



(Fh^h + 2Y He B Be ) 

2 {Y R B U + 2Y Rc B Rc ) MeV cm" 1 , 



(CI) 



where B B and B^ are of the form 

B = In 



7^(7-1) 



m e vc 



1 /v\ 2 

-J-) +K7,*). 



(C2) 



In the above equation, AE is the average excitation energy and AE = 15 and 41.5 eV for 
H and He atoms, respectively. In Eq. ( IClj) we have ignored the contributions from the He + 
ions as their abundance is much smaller than the abundances of H and He atoms. Note that 
the energy loss due to excitation and ionization of the bound electrons in atoms and ions is 



significant only when the positron kinetic energy is E^ 1 ^> AE. 

In the relativistic regime, three additional processes may be considered for positron energy 
loss. We follow the discussion in Ref. 18|] and first consider bremsstrahlung. In general the 
energy loss rate in an ionized plasma differs from that in neutral matter. However, Ref. [3] 
showed that the rates are the same for these two cases for 7 < 10 2 . The relevant energy loss 
rate for our problem is 



fdE e 



\ dx 



4a 



brem 



m e c A 



2.10 X 10~ n 7 



E e+ (2Y p + 6Y a 
ln(2 7 ) - - 



ln(2 7 ) - - 



MeV cm" 1 . 



(C3) 



Positrons can also lose energy through Compton scattering on the photons in the radiation 
field of the stellar surface region. The corresponding energy loss rate is 



(dE e 



V dx 



) = ^ f-^V) 2 a rad T 4 7 2 = 2.30 x lO^V MeV cm" 1 , (C4) 

/ Comp J \ m eC / 



where a ra d is the radiation constant. In the presence of a magnetic field B, positrons can 
lose energy through synchrotron radiation. This process is similar to Compton scattering 
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in that B can be viewed as a source of virtual photons. The relative importance of these 
two processes can be gauged by comparing the energy densities in the magnetic and the 
radiation fields: 

B2/m ~ 0.959 ( « V. (CB) 



a rad T 4 V100G, 
The energy loss rates — (dE e + /dx) eXtP \, — (dE e +/dx) CXt3i t, — (dE e + / dx)brem, and 
— (dE e + /dx)comp are compared in Fig. HI Note that the last two rates increase with in- 
creasing 7 (at least for the positron energies of interest here) whereas the first two rates 
decrease with increasing 7. Therefore, if the positron energy loss through bremsstrahlung 
and Compton scattering is unimportant in the relativistic regime, it can also be ignored in 
the nonrelativistic regime. Note also that — (dE e + /dx) eXjP \ exceeds —{dE e +/dx)c oxnp by a 
factor of at least ~2x 10 5 . So the positron energy loss through synchrotron radiation can 
be ignored for B<4x 10 4 G. We assume that B <C 4 x 10 4 G in the surface region of the 
star under consideration. 



APPENDIX D: CROSS SECTIONS FOR DIRECT ANNIHILATION OF 
POSITRONS AND POSITRONIUM FORMATION 

As an example of direct annihilation with bound electrons, we give the cross section 0"^°^ 



for annihilation of a slow positron with the electron in the H atom 

.2 \ 2 



°"da°,H = TT^cflf.H ( ] ~, (Dl) 

\m e c z J v 

where Z e ff,H is a function of the positron velocity v and can be approximated as 

Z e g h « 8.868 - 7.838 (—) - 102.77 (-Y + 527.38 (-Y 

vac/ vac/ \ac/ 

/ v \ / v \ 5 / v \ 6 

-978.68 1 — I + 773.15 1 — ) - 197.17 1 — ) . (D2) 
vac/ vac/ vac/ 



The above approximation of Z c s^n is valid for v < 0.7ac, which corresponds to positron 
kinetic energy of E^ < In/2 = 6.8 eV, i.e., below the threshold for Ps formation with 
the electron in the H atom. The cross sections crf°J, crf&ni anc ^ a daHe Q] are compared in 
Fig. [51 

The relative importance of direct annihilation and Ps formation depends on E^. Direct 
annihilation is dominant in the relativistic regime. As the positron becomes more and more 
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nonrelativistic, Ps formation becomes more and more important. The cross section for a 

nonrelativistic positron to form a Ps with a free electron at rest is (see e.g., floL [h]]) 

o7 / fc \ 2 00 / 
2 TTCt I Tb \ sr^ 9n/n 



CPs,f 



33/2 



h 



\ u(u+ 1)' 



(D3) 



where u = n 2 E^/Iji and g n is the Gaunt factor for forming the Ps in the energy state with 
pri ncip al quantum number n. The Gaunt factor is close to unity and can be approximated 



as 



20] 



g n ss 1 + 0.1728 



u 



1 



n 



2/3, 



u 



!)2/3 



0.0496 



u 2 + (4/3)u + 1 

n 4/3( u+1 )4/3 _ 



(D4) 



The calculation of the cross sections for Ps formation with bound electrons is rather complex. 
Here we use the measured cross sections <7p s h and <Tp s He for Ps formation with the electrons 

n n ' 

in the H [13J and He [14J atoms, respectively. The cross section for direct annihilation with 
free electrons multiplied by the number of free electrons per nucleon in the stellar surface 
region under consideration, Y e -af°J, is compared in Fig. [6]with the corresponding quantities 
^e-Cps,f) ^hCps,h, and VkefPs.He f° r Ps formation with free electrons and the electrons in the 
H and He atoms, respectively. 
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